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We investigate the effect of gravitational lensing by matter distribution in the universe on the cos- 
mic microwave background (CMB) polarization power spectra and temperature-polarization cross- 
correlation spectrum. As in the case of temperature spectrum gravitational lensing leads to smooth- 
ing of narrow features and enhancement of power on the damping tail of the power spectrum. 
Because acoustic peaks in polarization spectra are narrower than in the temperature spectrum the 
smoothing effect is significantly larger and can reach up to 10% for I < 1000 and even more above 
qq | that. A qualitatively new feature is the generation of B type polarization even when only E is 

■ intrinsically present, such as in the case of pure scalar perturbations. This may be directly observed 

^\ ' with Planck and other future small scale polarization experiments. The gravitational lensing effect 

is incorporated in the new version (2.4) of CMBFAST code. 
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I. INTRODUCTION 



Over the next few years a number of ground based, balloon and satellite experiments will measure CMB sky with an 
unprecedented accuracy and detail. The promise of a one percent precision on the measured power spectrum of CMB 
anisotropics requires a similar accuracy in the theoretical predictions, if we are to exploit all the information present 
in the data. The rewards will be rich: among other things this will allow an accurate determination of a number 
of cosmological parameters and testing of current structure formation theories |Q. In principle such a program is 
possible, since the anisotropics were produced when the universe was still in the linear regime, which makes the 
£f) • calculations of model predictions very accurate. In practice there are a number of important effects that need to be 
included if this goal is to be realized. One of the most important among these is the gravitational lensing effect. 

As photons propagate through the universe from their last scattering to our detectors they are randomly deflected 
by the gravitational force exerted upon them by the inhomogeneous mass distribution. Previous work has shown that 
• gravitational lensing has an effect on the temperature anisotropy power spectrum which is not insignificant . The 
Or random deflections smear out the sharp features in the correlation function or power spectrum, leading to a suppression 
Q ■ of acoustic oscillations. Gravitational lensing can also enhance power on the damping tail, causing it to decay less 
rapidly than predicted on very small angular scales || . Gravitational lensing effect on the temperature anisotropies 
& has been discussed several times in the literature and the formalism to calculate it using the evolution of density power 
spectrum both in linear and nonlinear regime has been presented in (^] . In this paper we extend this calculation to the 
two linear polarization power spectra and to the cross-correlation spectrum between temperature and polarization. 
Because acoustic oscillations are narrower for polarization spectra than for temperature, one expects gravitational 
lensing effect to be more significant in the former and indeed our results confirm this. In addition, a qualitatively new 
effect is the mixing between E and B types of polarization, which changes the pattern of polarization. The outline 
of the paper is the following. In §|| we develop the formalism: this section contains all the main analytic expressions 
needed for a numerical implementation of the effect. These have been numerically implemented in the new version of 



CMBFAST package (version 2.4) and require only a marginal increase in the CPU time for their evaluations. In § [II 



we compute the effect for a typical cosmological model and address the question of direct observability of the effect. 



We present the conclusions in §IV 



II. TWO-POINT CORRELATORS IN THE PRESENCE OF LENSING 



The large scale density fluctuations in the universe induce random deflections in the direction of the CMB photons 
as they propagate from the last scattering surface to us. This alters the power spectrum of both the temperature and 
polarization anisotropies. The quantity is responsible for the deflections is the projected surface density. Since the 
structures are not very correlated on large scales the gravitational lensing effect is only relevant at the small angular 
scales in the CMB. We may hence use the small scale limit formalism [0, which simplifies the calculations. 
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The observed CMB temperature in the direction is T(0) and equals the unobservable temperature at the last 
scattering surface T(0 + 60), where 68 is the angular excursion of the photon as it propagates from the last scattering 
surface until the present. In terms of Fourier components we have 

T{6) = f(0 + 66) 

= (2vr)- 2 J d 2 l e l[6+s6) T(l). (1) 

The same relation applies to the two Stokes parameters Q and U that describe linear polarization Q. In the small 
scale limit we take the direction of observation to be near z and we orient the local coordinate system orthogonal 
to z to define Q and U: Q is the difference between the photon intensities along x and y, while U is the difference 
between photon intensities along the two diagonals. In terms of Fourier components these can be expressed with 
opposite parity Fourier components E(l) and B(l) j/j 

Q(0) = Q{0 + 60) 

= (2tt)- 2 f d 2 l e d < 0+s °) [E(l) cos(2<^) - B(l) sin(2<^)] 



U{8) = U(0 + 60) 

= (2vr)- 2 J d 2 l e tl(e+se) [E(l) sin(2<^) + B(l) cos(2^)]. (2) 



The Fourier components satisfy 

{X{l)X{l')) = {2^f C kl 5 D {l-V), (3) 

with X = T, E, B and the average is over different realizations of the CMB field. 

The correlation function of the temperature between two points in the sky, and , only depends on their 
angular separation 9. With our choice of coordinate system the correlations between the polarization variables are 
also a function of azimuthal angle <f> fij . This is not the natural coordinate system in which to define these correlation 
functions. Instead one should align the local coordinate system with the great circle that connects the two points 
|| . With this choice parity conservation requires T and Q to be uncorrelated with U and the correlation functions 
only depend on separation 9. To obtain this set of correlation functions in our coordinate system we can calculate 
the correlation between the variables at the origin and another point separated by an angle 9 along the x axis. The 
correlation functions are calculated from equations ^ and (S), 

c T (e) = [ -^L ^wcos^ {ei i { so A -so B } } c 
J V'T^Y 

Cq{0) = Il^f e^ COB *(e^ A - 50B ) ) [C^cos 2 (2^) + C i3i sin 2 (2^)] 

Cu(e) = J ^ leC °^ l ^ liSGA - 5eB) ) [C £; sin 2 (2^) + C i|; cos 2 (2^)] 

Cc{9) = S e iiecoa ^{e il<sGA - 5e ^ ) C ei cos(24>i). (4) 

The remaining average in equation (||) is over the lensing fluctuations. Only the cross-correlation between Q and T 
is different from zero and we denote it Cc{&)- Even in the presence of lensing U does not cross correlate with either 
T or Q because these quantities have opposite parities. 

The correlation function of the excursion angle can be used to calculate the expectation value in equation (Q) ]3j , 

(exp^(<50 - 60')} ) = cxpj-^K 2 ^) + cos(2^)a 2 (fl)]} 

~l-^0) + cos(2^y 2 2 (0)]. (5) 

The exponential above has been expanded out assuming l 2 ct2{9) and l 2 cr 2 (0) are small and we numerically verified this 
to be an excellent approximation (see also |To| ) . The two functions characterizing the rms dispersion of the photons 
are |] 
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a 2 {0) = l&n 2 I k 3 dk r" P (fc,r = r o -x) W 2 ( Xl Xrec) [1 - J o (k0 sm K X )} 



(4(0) = 16tt 2 / k 3 dk I ^° P^(k,r = T - x) W 2 ( X ,Xrec) J 2 (k0 sm KX ). (6) 
Jo Jo 

We denote with X — T o ~ T the comoving radial distance, r is the conformal time and t corresponds to its value 
today. The comoving angular diameter distance is sin^- x — K~ x / 2 sini'T 1 / 2 ;^ in a closed universe [K > 0), x in a flat 
universe (K = 0) and (— K)~ x l 2 sinh(— K) x / 2 x, in an open universe (K < 0). The curvature K can be expressed using 
the present density curvature parameter fix = 1 — £l m — &a and the present Hubble parameter Hq as K = —Q,kH 2 . 
The power spectrum of the gravitational potential at time r is P{k,r) and W( X i Xrec) — ^ n K{Xrec — x)l sm if(Xrec), 
where Xrec is the radial distance to the last scattering surface at recombination. We include the non linear evolution 
of the power spectrum using the fitting formulae of Peacocks and Dodds || . These expressions are thus valid in a 
general Robertson- Walker metric, both in linear and nonlinear regime. The nonlinear effects are however not very 
important except on small angular scales. 

We use equations (||) and @ together with (|^) to obtain the final expression for the correlation functions, 

r idi i 2 i 2 

Ct(9) =J^C il {J (W)[1 - -a 2 {0)] + -„ 2 (9)J 2 {W)} 

r idi i 2 i 2 

C Q {9) + Cu(6) = J — (C El + C Sl ) {J (W)[1 - -ol{e)] + -a 2 2 (0)J 2 (l0)} 

r idi i 2 i 2 

c Q {6) - cu(0) = J — (c El - c Sl ) {uw)[i - ^im + joimMm + Mmn 

r idi i 2 i 2 

Cc{9) {J 2 {10)[1 - -a 2 (0)} + -a 2 (0)[J o (l0) + J 4 (10)}}. (7) 

The power spectrum in Fourier space has become the most widely used way to characterize the CMB anisotropies. 
We can compute the different power spectra from the correlation functions, 



Cti 


= 2tt 


Cei 


= 2?r 


Cbi 


= 2tt 


Cci 


= 2tt 



9d0 C c {0) J 2 {10)- (8) 
Equations (0) and (g) give the mapping between the observed CMB power spectra and the primordial one. Explicitly, 



C T i - Cft + W l n C TV 



Cei = C El + ^-[W l u + C El , + -[w( t - W l 2 \] C Bl , 



Cbi = C Bl + - C El , + \[W l u + W l 2l ] C Bl , 

C C i = C ei +W^C ei „ (9) 
the sum over I' is implicit. The window functions are defined to be 

Wfi = y J 9d0 J O (10) {a 2 (0)J 2 (l'0) - a 2 o (0)J o (l'0)} 

wli = l -j£ m Ji{w) {^ 2 2( ( >)lMi'e) + Mi'e)}~a 2 (0)j i (V0)} 

= \ [ 6de Mm {\a 2 (0)[J o {l'0) + J i {l'0)]-a 2 (0)J 2 (l'0)}. (10) 
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Equations @ and (0) are the main result of the paper. The results for polarization are new and represent a 
generalization of previous results for the temperature || . The important qualitatively new feature is that lensing 
mixes E and B polarization modes. On small scales where the lensing effect is important all cosmological models 
proposed so far predict only E type polarization and Cg l — 0. Lensing will however generate B type polarization 
in the observed field, Cbi = — W^] Cei' ^ 0. In the next section we will calculate the lensed power spectra 

in a typical cosmological model to address the significance of the effect. Before addressing this issue let us explore a 
simplified model to understand why the two polarization types arc mixed through lensing. 

So far we have introduced E and B type polarizations in Fourier space. We can also define real space quantities 



E{9) = (2tt)- 2 J d 2 l e M E(l) 

B{0) = (2tt)- 2 J d 2 l e M B(l). (11) 



These two quantities describe completely the polarization field and it proves easier to understand the effect of lensing 
in terms of these. They can also be computed directly from Q and U in real space fyl , 



E{9) = J d 2 9' uj{\e' - 0|) Q r {8') 

B(0) = J d 2 0' lu{\6' -6\) U r (6'). (12) 



We have defined Q r and U r , the Stokes parameters in the polar coordinate system centered at 0. If = then 
Q r = cos20' Q(0') - sin 20' U(0') and U r = cos20' U(9') + sin20' Q(6'). The window is w(6) = 1/6 2 (8 ^ 0), 
w{0) = (0 = 0). 

We consider as a toy model an unlcnscd polarization field which is a radial pattern around the origin in a ring of 
size O . We have Q(0) = E 9 5 D (6 - 8 )/2 cos 20 and £7(0) = E 6 6 D {9 - o )/2 sin 20. In the absence of lensing 
we would observe at the origin E(0 = 0) = Eq and B(9 = 0) = 0. This follows from U r = (because polarization 
is radial) and equation (p^). We now consider what happens if we shift the position of the photons by a random 
angle of size 58. Each segment of the ring will be mapped to a different position with angle (j> where \6<f>\ ~ \S9\/8q. 
Each segment acquires a random component of non radial polarization and since we are assuming the segment shifts 
are uncorrelated the integral of U r over <f> does not vanish. It has mean zero and variance (U 2 ) oc (59/8q) 2 Q 2 . The 
measured power in the B mode is (B 2 (9 = 0)) oc (S8/9o) 2 Eq. This example shows that the pattern of polarization 
vectors on the sky determines the amount of E and B type polarization. Lensing distorts this pattern by shifting the 
positions of the photons in the plane of the sky relative to the last scattering surface. It can thus generate B pattern 
out of initially pure E polarization. In the next section we discuss the amplitude of this effect. 



III. ESTIMATE OF THE LENSING EFFECT 



The effect of gravitational lensing on the temperature power spectrum has been studied in detail in previous work 
H . We will explore here the effects on polarization and asses the detectability of the B signal. In these examples 
we will use the cosmic concordance model E3| with Hq = 65 km s _1 Mpc _1 , Q\ = 0.65, Ct m = 0.35, f^/i 2 = 0.015 
and n = 1. The model is COBE normalized and satisfies the constraints on 8 h _1 Mpc scale from the abundance of 
clusters. 

Figure [l] shows the functions that characterize the rms deflection of the photons, <7o(0) and 02 (0). On small scales 
the rms relative deflection between two photons approaches 15%, justifying the weak lensing assumption over most of 
the sky (for caustic formation one requires 68 ~ 8, which will occur only on rare occasions). We also show ipA^iW^j] 
for I = 2000. The two window functions W{ z and are very similar, making their difference much smaller than their 
sum. One thus expects the generated B type polarization from a pure E type (and viceversa) to be very small. In 
fact Wg; is also very similar, the three windows only differ at the 1% level. The windows are oscillatory as a function 
of I' and their main contribution is concentrated around /. 

Figure ^| (upper panel) shows the lensed and unlensed power spectra. On intermediate scales lensing has the general 
effect of smearing the peaks in the spectrum by redistributing power. Polarization receives contribution only from 
velocity gradients at the last scattering surface and the acoustic peaks are very sharp. In contrast, temperature 
receives contributions both from velocity and density of photon-baryon plasma and the two are out of phase with 
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FIG. 1. The upper panel shows the two functions ao{9) and <T2(9) for the cosmic concordance model discussed in the text. 
The lower panels show ^[W^ ± W^] for I = 2000. The Window function for the mixing between E and B is is not well localized 
in I space and is rapidly oscillating as a function of I. 
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each other, leading to a partial cancellation of acoustic peaks. One thus expects the effect of lensing to be larger for 
polarization and indeed we find that it is approximately a factor of two larger for E polarization (figure 0, bottom 
left panel). In the damping tail the power is enhanced over the unlensed case as can be seen in the general trend in 
the relative difference between lensed and unlensed spectra (figure |^) . 

To asses the detectability of the induced B polarization we will focus on the Plack mission. We assume that 
measurements of the temperature, E type polarization and T — E cross correlation have allowed the determination 
of the cosmological parameters accurately enough so that we know the approximate shape of the B power spectra 
induced by lensing. We will attempt to determine just one parameter, the overall normalization of this B signal, 
which expresses the overall effect of the matter fluctuations along the line of sight. If the signal to noise were large 
enough we could also attempt to determine the power spectrum of the lensing effect, by exploring the signal in B as 
a function of I. 

The relative error on the overall amplitude of the induced B component /3 is 



A^ = / 2 

V /»*» £/( 2Z + !)/(! + w p l B- 2 /C B i) 2 



(13) 



f s ky is the fraction of observed sky which we take to be f s k y — 0.8. We assume the noise and beam width to be Wp 1 = 
(0.025 (J,K) 2 and B^ 2 = e' 2<T t with cr b = 9f whm /2\/2 In 2 = 9 x 10~ 4 . The effective number of modes contributing 
information is N eff = fskyJ^iC 21 + + Wp 1 B~ 2 /C B i) 2 , so that A/3//5 = y/2/N e /f. For the experimental 

specifications above N e ff w 8 and A/3//3 w 0.5. This means that there are only 8 independent modes of B type 
polarization that can be observed with Planck, each a gaussian random variable with mean and so the signal will 
be at the limit of detectability by Planck. The B signal peaks at I = 1000 and a ground based experiment observing 
a small patch of the sky for sufficiently long time to reduce the noise per pixel would be more effective to detect this 
signal. Gravitational lensing induced B is also not a significant contaminant of the B polarization expected from 
tensor modes in inflationary models the latter being dominant on large angular scales (figure ||, bottom right 
panel, where we assumed tensors and scalars are of equal amplitude in the temperature on large scales). 



IV. CONLUSIONS 



The change in the CMB spectra induced by random deflections of the photons by the large scale structure of the 
universe has to be included in the calculations of the anisotropies when comparing theory and observations. We 
have presented a formalism to calculate the effect of gravitational lensing on all four CMB power spectra. The final 
expressions are given in a compact and numerically efficient form that is adequate for numerical implementation. These 
expressions have been implemented in the CMBFAST package, freely available from the authors. The computational 
time does not significantly increase with this feature and since the effect can be significant we feel the lensing effect 
should be included in any calculation where high precision accuracy is important, such as in the design and analysis 
of the CMB experiments. The method is self-consistent in the sense that for any cosmological model we use the 
actual power spectrum computed in the code to compute the lensing effect. The power spectrum is normalized to 
COBE using the CMB spectrum computed from the same code output. This approach gives the correct amplitude 
of gravitational lensing effect for the particular model in question. It should be pointed out that if the model is 
inconsistent with the small scale constraints such as as normalization then the amplitude of the lensing effect will 
also be incorrect. For example, COBE normalized standard CDM has a much larger gravitational lensing effect than 
what we find using the cosmic concordance model, but this is only a reflection of standard CDM model having too 
much small scale power to be compatible with small scale constraints. For models that are correctly normalized on 
small scales the relative change of the polarization power spectra can reach 10% at I ~ 1000 and even more at higher 
I. This is larger than in the temperature spectrum and is caused by the sharper acoustic peaks in the polarization 
spectra. In the damping region the lensed spectra show an enhancement above the unlensed spectra just like in the 
temperature case, although the sensitivity to this effect is too small for satellite missions to detect it in polarization 
& 

Gravitational lensing also mixes E and B type polarization by deforming the polarization pattern on the sky relative 
to that at the last scattering surface. This will generate B type polarization out of the E type polarization even if 
there was no B present at the last scattering surface. This induced B mode is rather small in typical models and will 
be only marginally detectable by the Planck Surveyor. It peaks at fairly small angular scales around I ~ 1000 and so 
does not affect the measurement of gravity waves from B polarization on larger scales. A ground based experiment 
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FIG. 2. The upper panel shows the T, E and C power spectra. Dashed (solid) lines correspond to the lensed (un- 
lensed) spectra. The bottom left panel shows the relative difference between the lensed and unlensed spectra for T and E 
(SCi/Ci = (Cl ensed - c? nlensed )/Cr lensed ) and shows both suppression of oscillations and enhancement of power on small 
scales. The bottom right panel shows the B type polarization induced by lensing. We include the E and B spectra for 
inflationary model where scalars and tensors produce equal amount of power in the temperature on COBE scales (T/S = 1). 
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observing a small patch of the sky would be more suitable to observe this effect and would allow one to determine 
the power spectrum of matter fluctuations. This would allow a detection of a combined gravitational effect from 
structures spanning a much larger range in redshift space than currently reachable by other methods. 
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